Competition between fluctuations and disorder in frustrated magnets 
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We investigate the effects of impurities on the nature of the phase transition in frustrated magnets, 
in d = 4 — e dimensions. For sufficiently small values of the number of spin components, we find no 
physically relevant stable fixed point in the deep perturbative region (e <C 1), contrarily to what is 
to be expected on very general grounds. This signals the onset of important physical effects. 
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The influence of disorder in solid state physics is a very 
important problem. In many cases, it is expected that 
inhomogeneities, such as impurities, defects in the lattice 
structure, etc. may induce a completely different behav- 
ior in the system. Let us discuss two such situations in 
the context of phase transitions, which we shall be con- 
cerned with in this report. In systems undergoing second 
order phase transitions, disorder most often changes a di- 
vergence in the specific heat at the critical temperature 
into a cusp. The archetype of such a situation is the 
d = 3 Ising model, for which both experimental and the- 
oretical works show that the critical exponent a changes 
from positive to negative when disorder is added (for a 
review, see Ref. Q). This is well understood, thanks to 
the Harris criterion^ which states that for a large class 
of disordered systems^ the critical exponent v& - which 
governs the singularity of the correlation length - must 
satisfy the inequality v& > 2/d, where d is the number 
of spatial dimensions. Using the hyperscaling relation 
a = 2 — vd, this inequality translates into a < 0, which 
in turn implies that there is no divergence in the spe- 
cific heat. More generally, Harris criterion leads to the 
conclusion that adding disorder to a pure system with 
v p < 2/d has a dramatic effect, since it must change the 
critical exponents, and therefore the universality class of 
the phase transition. 

The influence of disorder on first-order phase transi- 
tions has also been considered for a long time. It was soon 
proposed that impurities tend to reduce the discontinu- 
ities of the first derivatives of the free energy (jump of the 
magnetization, latent heat). It was argued that these dis- 
continuities could even disappear in some cases^ Again 
disorder changes qualitatively the properties of the sys- 
tem, by inducing a second-order phase transition. Later 
on, this "rounding effect" of disorder was proven to take 
place in a large class of systems by Aizenman and Wehr 5 
(referred to as AW in the following). Up to now there is 
no experimental realization of this phenomenon. 

It was recently proposed by one of us& that such 
a second-order phase transition induced by impurities 
could actually be studied by experimental means in frus- 
trated magnets such as stacked triangular antiferromag- 
nets (CsMnBr3, CsMnl3, etc.) and rare-earth helimag- 
nets (Ho, Dy, Tb). In the pure case, these materials 
have been very much studied both experimentally and 
theoretically since the 1970s (see Refs. Q and |g for re- 



views). The early calculations indicate a second-order 
phase transition (see also Ref.0), in accordance with the 
observed power-law behavior of thermodynamic quanti- 
ties. However, there is a wide dispersion of the measured 
critical exponents. This lack of universality has been re- 
cently interpreted as a signature of a weakly first-order 
phase transitioniifliii This issue is still very debated. In 
any case, we believe that it would be very instructive 
to study the influence of disorder in such systems. On 
the one hand, if the transition is of second order in the 
pure case, one expects a dramatic change of the critical 
exponents when impurities are added. Following Harris 
criterion, the exponent v, which is always found to be 
smaller than 0.57, should jump above 2/3. This varia- 
tion is at least three times bigger than in the case of Ising 
magnets and might be easier to observe experimentally. 
On the other hand, if the transition is weakly of first 
order, it is expected to be turned into a continuous one 
in the presence of impurities, on the basis of AW. As a 
consequence, one would observe universal critical expo- 
nents, in contrast with the scattered values reported in 
the pure case. This would provide the first experimen- 
tal realization of the rounding effect of disorder. This 
last scenario is supported by a recent nonperturbative 
renormalization-group calculation in d = 3i* Actually, 
the results of AW are not restricted to d = 3, and apply 
to any dimension lower than four^i In the present re- 
port, we study the rounding effect of disorder within the 
standard e expansion. As we shall see, the perturbative 
results do not coincide with the conclusions of AW, sig- 
naling the onset of physical effects, yet to be determined. 

In the following, we describe the one-loop result, 
from which we argue that higher-order contributions are 
needed in order to get the qualitative behavior of the 
system. We show that the e expansion of the fixed-point 
coupling constants is singular for some physically relevant 
values of the number of spin components. Using a simple 
geometrical picture, we explain the origin of such singu- 
larities, and recast the perturbative series into a regular 
form. Our discussion is quite general, and can be applied 
to a wide class of theories displaying the same qualita- 
tive behavior. We then perform a two-loop calculation 
and analyze the effect of disorder using this method. 

Let us first describe the Hamiltonian which is used 
in the following. The theoretical studies of the phase 
transition in the materials considered here started in the 
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1970sfi2ii^i It was then realized that the relevant order 
parameter is not a simple vector (as for a ferromagnetic- 
paramagnetic phase transition) , but a set of two vectors 
(pi and tf2- This is a direct consequence of frustration, 
which makes the ground state more degenerate than in 
the corresponding nonfrustrated system. In physically 
relevant situations, (p\ and (p2 are taken to have two or 
three components, depending on the anisotropies of the 
crystal. Here we consider a direct generalization where 
the order parameter is composed of rn vectors, each hav- 
ing n spin components, and we gather these into a n x m 
matrix: <j) — (tpi, • • • , <p m )- Considering only the relevant 
"y 4 -like" terms, the hamiltonian reads: 



S{<f>) = I d d x[-Tr{d% 
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where p = Tr('00) and r = Tr( (fxf>tf>(f)). Adding impu- 
rities induces a spatial dependence of the coupling con- 
stants, through the local density of impurities. Using 
the replica trick, it is, however, possible to average over 
disorder, and to study the properties of the system by 
considering a translational invariant Hamiltonian. The 
price to pay is that we now have to deal with o replicas 
of the field, and to take the limit o — > at the end of the 
calculation. The < / 9 4 -like hamiltonian then reads^ 



i=i J i=i 



(2) 



where I indices the different copies of the field. Note that 
only the region u% > 0, mu\ + 2u2 > 0, 113 < is of 
physical relevance ii^ 

Let us describe the situation at the one-loop level. The 
(3 functions were derived in Ref. 0, and can be read off 
on Eq. below. Out of the eight solutions to the set of 
three equations /3j = 0, only four are of interest here. The 
situation of pure frustrated magnets is recovered when 
«3 = 0. In this plane, we retrieve the two associated 
fixed points, noted P + and P~ in the following. In ad- 
dition to these well-known fixed points, we find two new 
ones which are of particular interest in our study, and 
which we call D + and D~ . Let us discuss the behavior 
of the fixed point when n is varied, begining with P + and 
P~ H For large values of n, they have real coordinates 
in coupling-constant space. When n is decreased, they 
approach each other, and coincide at the critical value 
ri+ = 4(3 + V6) — 21.8 (throughout this report, the 
numerical values are given for the physically interesting 
case 171 — 2). Below rip , the two fixed points get com- 
plex coordinates and are not physical anymore. Finally, 
when n reaches n~ = 4(3 — \/6) — 2.2, P + and P~ 
reappear, however, not in a region of physical interest. 
These features follow directly from the expression of the 
coupling constants at P^, which read (i = 1, 2, 3) 



e Ai (n) ±eB l (n) \Jn-n pfi \Jn- 



(3) 



We observe the very same qualitative behavior for the 
couple Z? ± , and we denote by n^ Q the associated critical 



values of n. Surprisingly, the two couples of fixed points 
P and D annihilate for the same values of n (namely, 
rip q = rijQ = nj), and at the same location in coupling- 
constant space. In the following, we focus our attention 
on the vicinity of the upper critical value 71q , the only 
physically interesting one, and we drop the superscript. 

In summary, at one loop, when n > no, the fixed point 
P + is stable, and the disorder is irrelevant. For n < no, 
there is no stable fixed point, and the phase transition is 
likely to be of first order. This is very astonishing, being 
difficult to reconcile with the general results of AW. Thus 
one can expect that higher-order corrections will resolve 
this mismatch. In fact, we expect a qualitative change al- 
ready at two loops. Indeed, for the pure frustrated mag- 
nets, at the value of n where the two fixed points P 
meet, one has^ v = 2/d~e 2 (V6- l)/200 + o(e 3 ) < 2/d. 
Using Harris criterion, we conclude that the disorder 
should be relevant, at least for this value of n. This in 
turn suggests that the degeneracy observed at one loop 
should be resolved at higher orders, namely, that the two 
couples of fixed points P ± and D ± should annihilate in- 
dependently and for different values of n. 

Before embarking in higher-loop calculations, it is use- 
ful to do a short digression and to consider the phe- 
nomenon of annihilation of fixed points on more general 
grounds. For the sake of simplicity, consider a theory 
with two coupling constants u\ and U2- The points where 
the associated (3 functions vanish form two curves in the 
plane (u\,U2), whose intersections are the fixed points. 
When a parameter is varied - in our case the number of 
spin components - these curves move, and so do their 
intersections. A typical situation is represented in Fig. 
^ where one sees the annihilation of two fixed points for 
n = n c . This qualitative picture is already observed at 




(a) n > n c 



(b) n = n c 



(c) n < n c 



FIG. 1: Annihilation of two fixed points at n = n c . The 
lines represent zeros of the f3 functions and their intersections 
represent the fixed points. 



one loop and it is clear that it does not change when 
higher-order corrections are included: the curves of Fig. 
^are only slightly displaced. As a consequence, the posi- 
tion (u\, U2) where the two fixed-points meet, as well as 
the corresponding value n c of the number of spin com- 
ponents, receive small corrections at each order of per- 
turbation theory and therefore depend on e. From these 
considerations, it is an easy matter to deduce the generic 
form of the fixed-point coordinates near n c (e). Expand- 
ing the (3 functions in coupling constants up to second 
order around u*, one is left with quadratic expressions 
as in the one-loop case. The coordinates of the fixed 
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points then read 

uf = e Ri (n, e) ± e Si (n, e) y/n- n c (e) , 



(4) 



where R4 and Si are polynomials in e. The previous 
formula is very appealing when compared with the one- 
loop form ® in the vicinity of Uq , the main difference 18 
being the e dependence of the critical value of n . The 
whole previous argument is very general, and only makes 
use of the geometrical picture of Fig. ^ I n particular, 
it does not rely on perturbation theory and Eq. |QJ is, 
in this sense, a nonperturbative result. In the following, 
we make extensive use of Eq. Q which, when combined 
with perturbation theory, provides a very powerful tool. 

Let us finally make contact with the usual e expansion 
of the fixed-point coupling constants, the coefficients of 
which are found to be singular in the limit n —> rio^ 
We can easily understand the origin of such singularities 
in the light of our previous discussion: writing n c (e) = 
n c ,o + eiic.i + • • • and expanding Eq. Q in powers of 
e, one gets terms of the form e^/(n — no)" 1-1 / 2 , with 
m < £ — 1. Therefore, the regular form Q corresponds 
to a resummation of infinitely many singular terms. In 
practice, one reconstructs the polynomials Ri and Si by 
matching the e expansion of Eq. Q for n ^ no with the 
corresponding perturbative series. 

Comming back to the problem at hand, we now dis- 
cuss the expected influence of higher-loop corrections on 
the one-loop flow diagram. The critical values of n in- 
troduced before, now depend on e, and we shall refer to 
them as n p (e) and rid(e)- As already explained, Harris 
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FIG. 2: Flow diagram above and below n s . The fixed point 
D + crosses the plane U3 = through P + , and gets stable 
afterward. 

criterion indicates that disorder should be relevant, at 
least for n — n p (e). This can only happen if, for some 
particular value n s (e), either D + or D~ crosses the plane 
«3 = through P + , and exchanges its stability with it. 
Moreover, following AW, we expect that for n < n p (e) 
the phase transition is governed by a disordered fixed 
point. This favors the first scenario, which is depicted in 
Fig. |21 In addition to n p (e), rid(e), and n s (e), we shall be 
interested in a fourth particular value of n, noted nn(e), 
for which the critical exponent v associated with P + just 
equals 2/d. Exploiting the Harris criterion and the topol- 
ogy of the flow, we see that n s is larger than n p , rid, 
and Tin- 

With this in mind, we perform a two-loop calculation 
of the P functions. Using dimensional regularization and 



minimal subtraction, we get 

fh = -uxe + 2u 4 {ui 2 / 8 ^,i + 4wi(3u 3 + u 2 hs,o) 
+ 12m 2 2 } -W{3u 1 3 I 14fi . 3 + 2u 1 2 (22u 2 I 1 . lfi 

(5a) 

+«3/58,0,ll) + Ml(68u 2 U3^1,l,0 + 2li 2 47,5,5 

+w 3 2 J82,5) + 24u 2 2 (6u 3 + 142/3,1,0)}, 

02 = -u 2 e + 4u 2 w 4 {(6ui + U2-?4,i,o + 6u 3 } 
-4u 2 w 4 2 {6m 2 2 /i7, 3 ,i + 4u 2 (ui + U 3 )/29,u,o (5b) 

+U\{UI + 2u 3 )/ 82:0 ,5 + M3 2 ^82,5}, 

03 = ~u 3 e + 2u 3 v 4 {2uil 2 fl,i + 4u 2 /i,i,o 

+■"3^8,1} ~ 4m 3 u 4 2 {10m 2 2 / 3 4,i -I- 3m 3 2 Ji4, 3 (5c) 
+4m 2 (5mi + 111*3)11,1,0 + ui(5«i + 22m 3 )/ 2j0 ,i}, 

where w 4 = 1/ (327r 2 ), Iij.k = i + j( m + n ) + k(mn), and 
Jij = i + j(mno). We also give the expressions of the 
critical exponents: 



(6a) 



77 = 2u|{l 2 ,o,iUi(ui + 2u 3 ) + 4/14,0^2 (ui + u 3 ) 

+ 2I 3AA u 2 2 + 2J 2A u 2 3 } , 
l/v = 2- 2v 4 {/ 2 ,q,iUi + 2/i,i, M2 + ^2,1^3} + 577 (6b) 



By using Eq. we can safely study our expressions 
near n = hq. The first surprising result is that the four 
values of n introduced above are still degenerate at order 
e, i.e. u c ,i = ni = -2(18 + 7%/6)/3 ~ -23.4, for c = 
d, p, s, H. In principle, the next correction, i.e. n. c ,2, 
is obtained by a three-loop calculation. However, due 
to Eq. Q we can show that the degeneracy is resolved 
at order e 2 without performing explicitly a three-loop 
calculation. Let us illustrate this point on a simple case. 
When n — n s (e), the fixed points P + and D + coincide 
and the coupling constant u 3 of D + vanishes. Inserting 
n = n s (e) = uq + eri\ + e 2 n. s ,2 in Eq. Q, and expanding 
the equation u\ = to second order in e, one obtains 



n d , 2 



(d e + nid n )R 3 (n, e)|n=n , e =o 
S3 (no, 0) 



(7) 



The expressions appearing in the above equation involve 
only two-loop expressions and can therefore be computed 
here. For m = 2, we find 



™ S ,2 - n d ,2 



(1092 + 463\/6)/50 ~ 44.5. 



(8) 



Using the same technique for the equations Ui(D + ) = 
Ui(P + ) at n = n s (e) (i = 1,2), as well as v{P + ) = 2/d 
at n = fi«(e), we get 

n B , 2 - n p , 2 = n H , 2 - ^ P ,2 = (7V6 - 12)/50 ~ 0.103 . (9) 

We see that n s (e) = nn(e) up to second order: the Har- 
ris criterion is satisfied precisely at this value of n where 
disorder becomes relevant. We mention that a three-loop 
calculation for the pure (1*3 = 0) frustrated magnets 14 
yielded the value: n Pj2 ~ 7.09. Combining this with our 



4 



"two- loop" results © and ©, we are able to compute 
the "three-loop" contributions: n^ 2 = —37.3, n s 2 — 
bh.2 = 7.19. We find the same features for any value 
of to, except for the numerical results. For example, for 
?7i = 3, we get 7id,2 = —48.2 and n Sj 2 = nu,2 = 11.1. 
We summarize our results in Fig. There are essen- 
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FIG. 3: Different types of phase transition in the (n, d) plane. 
The lines represent n s (e), n p (e) and n<j(e) from top to bottom. 
The splitting between n s (e), and n p (e) has been magnified 

tially three distinct regions in the (n, d) plane. In the 
first one [n > n s (e)), the transition is of second order 
and is governed by the pure fixed point P + . In the sec- 
ond one [ rid ( e ) < n < n s( e )L the disorder is relevant 
and the system undergoes a second-order phase transi- 
tion governed by the fixed point D + . This region is di- 
vided by the curve n p (e) (dotted line in Fig. |3J. Above 
it, P + still exists, but is unstable against disorder, as ex- 
pected from the Harris criterion. Below, the pure fixed 
point has disappeared: the fluctuation-driven first-order 
phase transition observed in the pure system is turned 
into a continuous one by impurities, in line with AW. 
In the whole second region, the disorder induces a new 
universality class, and therefore new critical exponents. 
For instance, on the line ?7.d(e), one has 77 = e 2 /48 and 



v = 2/d + e 2 (13 + 7\/6)/200 up to o(e 3 ). We note that 
the experimentally relevant cases {d = 3, n = 2,3) be- 
long to this region. Although our results can only serve 
as a rough indication in d — 3, they are corroborated 
by the conclusions of a nonperturbative approach, and 
support the recent proposal to study experimentally the 
rounding effect of disorder in frustrated magnets £ Fi- 
nally, there is a third region [n < nd(e)], where we find 
no physically relevant stable fixed point: here, the dis- 
order does not counterbalance the effect of fluctuations, 
and the phase transition most probably remains of first 
order. We stress that for n < no, the deep perturba- 
tive regime e -C 1, where our approach is under control, 
belongs to this third region. 

If the results found in the first two regions are well 
understood in terms of the Harris criterion and AW, 
the situation is still puzzling in the third region. In- 
deed, perturbation theory predicts a discontinuous phase 
transition where AW state that the latent heat vanishes, 
which in turn strongly indicates a continuous transition. 
We suggest three possible scenarios which could explain 
this mismatch. It could happen that, although the la- 
tent heat vanishes, the phase transition is still of first 
order—. However, simple Landau-type arguments indi- 
cate that this is not very likely. Another possibility is 
that the replica symmetry, which was implicitly assumed 
to hold, may actually be broken. This could be inves- 
tigated following the lines of Ref. Finally, the ex- 
pected second-order phase transition could be governed 
by a nonperturbative fixed point, unreachable by the ap- 
proach we considered here. These three possibilities are 
very interesting and deserve more investigations. 
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Strictly speaking, this holds if the number of spin compo- 
nents is larger than two 

We mention that in principle one can repeat this construc- 
tion near n = , so as to get a form with two square 
roots, similar to J!|J. 

There could be a jump in the magnetization, as observed 
in Ref.El 



